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Abstract: We present a detailed study of a specific class of graph that can potentially 
contribute to the proton-proton double parton scattering (DPS) cross section. These are 
the '2vl' or 'single perturbative splitting' graphs, in which two 'nonperturbatively gen- 
I erated' ladders interact with two ladders that have been generated via a perturbative 

■ 1 — )■ 2 branching process. Using a detailed calculation, we confirm the result written down 

D , originally by Ryskin and Snigirev - namely, that the 2vl graphs in which the two nonper- 



turbatively generated ladders do not interact with one another do contribute to the leading 
order proton-proton DPS cross section, albeit with a different geometrical prefactor to the 
one that applies to the '2v2'/'zero perturbative splitting' graphs. We then show that 2vl 
graphs in which the 'nonperturbatively generated' ladders exchange partons with one an- 
other also contribute to the leading order proton-proton DPS cross section, provided that 
this 'crosstalk' occurs at a lower scale than the 1—7-2 branching on the other side of the 
graph. Using results from past studies of crosstalk interactions in the context of shadowing 
corrections to DIS, we argue that even though crosstalk interactions are suppressed by 
colour effects, the overall impact of crosstalk effects on the DPS cross section may be sig- 
nificant. We give an explicit formula for the contribution from the 2vl graphs to the DPS 
cross section, and combine this with a suggestion that we made in a previous publication, 
that the 'double perturbative splitting'/'lvl' graphs should be completely removed from 
the DPS cross section, to obtain a formula for the DPS cross section. It is pointed out that 
there are two potentially concerning features in this equation, that might indicate that our 
prescription for handling the Ivl graphs is not quite correct. 
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1 Introduction 

We define double parton scattering (DPS) as the process in which two pairs of partons 
participate in hard interactions in a single proton-proton (p-p) collision. This process is 
formally suppressed relative to the usual single parton scattering (SPS) mechanism by 
A^/Q^, where is the scale of the hard scattering(s) involved, and is the scale of 
nonperturbative physics. However, DPS can contribute important backgrounds to SPS 
processes that are suppressed by small or multiple coupling constants, such as Higgs or 
new physics signals [1-4]. It is also an interesting process to study in its own right, as it 
reveals novel information concerning correlations between partons in the proton. For any 
final state AB which can potentially have been produced via two independent scatterings 
yielding A and i?, there is a region of final state phase space in which the DPS production 
mechanism is competitive with the SPS one - namely, the region in which the transverse 
momenta of A and B are small [5, 6]. This fact offers some scope to study DPS in detail - 
indeed it has been used in past experimental extractions of DPS [7-12] . Finally, the rate 
of DPS relative to SPS for a given final state AB increases with collider energy, as lower 
X values are probed where the population of partons is larger. This means that DPS will 
be more important at the LHC than at any previous collider. For these reasons, there has 
been a considerable increase in interest in the phenomenon of DPS in recent years from 
both the experimental and theoretical communities, and four international workshops have 
been organised on the theme [13-17]. 

If one assumes that the two hard processes A and B may be factorised, then the total 
cross section for the process — )■ AB + X via double parton scattering should be of the 
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following form: 



°^ X] 7^i:^]^dxaTij{xi,X2,r;Q\,Q%)rki{x3,X4,-r]Q\,Q%) (1.1) 

i,j,k,r ^ a=l 

xaik^A^s = xiX3s)aji^B{s = X2X4S) 

oc ^ / (fbY\_dxaTij{xi,X2,b;Q%Q%)rki{x3,Xi,b;Q\,Q%) 

i,j,k,l a=l 

xaik-,A{s = xiX3s)aji^Bis = X2Xis) 

The a symbols represent parton-level cross sections. ry(xi, X2, &; Q^, Q^) is the impact- 
parameter space two-parton GPD (b-space 2pGPD), whilst ^ijixi, X2,r;Q\,Q'^) is the 
transverse momentum space 2pGPD (r-space 2pGPD). Tij{xi,X2, b; Q\, Q^) has a proba- 
bility interpretation as the probability to find a pair of quarks in the proton with flavours 
ij, momentum fractions X1X2, and separated by impact parameter b, at scales Qa and Qb 
respectively [5, 6, 18]. The r-space 2pGPD is the Fourier transform of this with respect 
to b, and has no probability interpretation, r is related to the transverse momentum im- 
balance of one of the partons emerging from the proton between amplitude and conjugate 
(for more detail, see Section 2 of [5]). 

Since the experimental extraction of DPS relies on the fact that the DPS cross section 
differential in the transverse momenta of A and B, qa and qs-, is strongly peaked at small 
Qa and qB, it is perhaps the DPS cross section differential in qA and qs rather than 
the total cross section that is more relevant for making experimentally testable predictions 
[5, 6, 18, 19]. In the region of qA,qB of interest (i.e. q\,q\ ^ Q\-,Q\) this quantity is 
described in terms of transverse momentum dependent 2pGPDs (TMD 2pGPDs), rather 
than the collinear 2pGPDs appearing in (1.1). On the other hand, it is expected that for 
h? <^ q\^q\ <^ Q\iQ\ the TMD 2pGPD should be expressible in terms of the collinear 
2pGPD and a perturbatively calculable piece [5, 6]. In that case there is a 'collinear part' 
of the differential cross section whose structure closely resembles the total cross section 
formula (1.1). Knowledge of how the total DPS cross section is to be treated should be 
helpful in establishing the correct way to treat this collinear part. It is with this ultimate 
purpose in mind that we continue to discuss only the total cross section for DPS in the 
remainder of the paper. 

There are several classes of graph that can potentially contribute to the leading order 
(LO) p-p DPS cross section (we restrict our attention to leading order, or leading loga- 
rithmic accuracy, in this paper). These are sketched in figure 1. Note that the partons 
emerging from the grey proton blobs in the figure are nonperturbatively generated partons 
- i.e. ones existing at a low scale ~ A.qcd - and that we've taken all the hard processes in 
the figure to be the production of an electroweak gauge boson with positive invariant mass 
(denoted by a wiggly line). We'll refer to the different classes of graph (a), (b), and (c) as 
2v2, 2vl and Ivl graphs respectively, for obvious reasons. 

In the paper [20], we carefully examined graphs of the Ivl type. We found that the 
treatment of these graphs by a long-established framework for calculating the p-p DPS 
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Figure 1. Some types of graph that can potentiahy contribute to the DPS cross section. 

cross section [21-23] was unsatisfactory, and suggested that no part of these graphs should 
be included as part of the leading order p-p DPS cross section. Since then, this suggestion 
has also been made in a number of other papers [19, 24]. 

In light of this discovery, a careful re-analysis of other classes of graph that can poten- 
tially contribute to the LO DPS cross section would seem appropriate. In this paper we 
will pay particular attention to the 2vl graphs in which there is only a single perturbative 
splitting, such as that drawn in figure 1(b) (we'll also discuss to a certain extent 2v2 graphs 
such as 1(a) in which there are no perturbative splittings, although it should be reasonably 
clear that these should be included in the LO DPS cross section). 

In section 2, we will begin to address the issue of whether contributions from the 2vl 
graphs should be included in the LO DPS cross section, and what form these contributions 
should take. We'll do this using a similar strategy as we employed for the Ivl graphs in 
the paper [20]. That is, we'll take a 2vl graph with the simplest possible structure (i.e. 
the structure of figure 2) and see whether there is a 'natural' part of the cross section 
expression for it that is proportional to 1/Rp {Rp = proton radius), and also contains 
a large logarithm associated with the 1—7-2 splitting. The large logarithm should be 
associated with transverse momenta of the partons emerging from the 1 — )• 2 splitting 
being <^ (where we take Q\ = = for simplicity). If there is such a structure in 
the 2vl graph, then this part of this graph should be included in the LO DPS cross section. 
Furthermore, if there is a log((5^/A^)/i?p structure in the simplest 2vl diagram, then we 
expect there to be a log((5^/A^)"/iip piece in the more general 2vl diagram of figure 1(b) 
that should also be included in the LO DPS cross section. This will be associated with 
the branchings in the diagram being strongly ordered in transverse momentum. From the 
structure of the contribution to the LO DPS cross section coming from the simplest 2vl 
diagram, we'll be able to write down a resummed expression for the contribution to the 
LO DPS cross section coming from 2vl diagrams with the structure of figure 1(b). 

The results that we obtain in section 2 have in fact already been written down in 
the papers [19, 24-26], and one can view the content of that section as a more detailed 
re-derivation of some of the results in those papers. In section 3, we will however estab- 
lish a further result with regard to the contribution of 2vl graphs to the LO DPS cross 
section. We will discover that the final formula that we obtained in section 2 is incom- 
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Figure 2. The simplest structure possible for the 2vl graph. 



'Crosstalk interactions' 
between ladders 




Usual' ladder rungs 



1 — > 2 branching scale, k 



Figure 3. Generic 2vl diagram including 'crosstalk' that we argue contributes to the 2vl DPS 
cross section at the leading logarithmic level. 



plete, and that there are further diagrams of the 2vl type that contribute to the LO DPS 
cross section. These diagrams involve non-diagonal crosstalk interactions between the two 
nonperturbatively generated parton ladders, at scales lower than the perturbative 1 — )• 2 
ladder branching on the other side - an example diagram of this type is sketched in figure 
3 (note that there are no such crosstalk interactions in figure 1(b)). This result is again 
estabhshed by analysis of the simplest Feynman graphs of the appropriate type. In section 
3 we'll also discuss in detail the issue of colour in relation to the crosstalk interactions. 

In section 4, we combine the results from sections 2 and 3 with our suggestion from 
[20] with regards to the Ivl graphs, to give a suggested expression for the LO cross section 
for DPS. We do however point out two potentially concerning features in this equation 
that may indicate that the suggestion of [20] to completely remove the Ivl graphs from 
the DPS cross section is not quite the correct prescription. 



2 'Two versus One' Contributions to the DPS Cross Section 

In this section, we show explicitly that for a 2vl diagram with the structure of figure 2, 
there is a part of the cross section expression that contains a DGLAP-type large logarithm 
and a factor of order 1/i??, which should be considered as part of the LO DPS cross section. 
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Figure 4. (a) An example of a '2vl' DPS-type scattering diagram, (b) An example of a '2v2' 
DPS-type scattering diagram. The thick grey lines are protons, whilst the grey circles are proton 
vertices. The labels on the lines correspond to the four momenta of those lines. 

We present details of the calculation only for the particular flavour-diagonal contribution 
to the gp — )• gqq + X — )• 7*7* + X process presented in figure 4(a), where the two off-shell 
photons both have a positive invariant mass. However, the general method outlined below 
can be applied to any diagram of the appropriate structure, and will always give a large 
logarithm provided that the corresponding process is allowed in the collinear limit (apart 
from issues of Jz nonconservation at the splitting vertex). 

In the calculation of the cross section for figure 4(a), we will have to include a wave- 
function factor or hadronic amplitude (p to find two nonperturbatively generated partons 
in the proton, at the amplitude level in the calculation. It would be inappropriate to try 
and calculate a 2vl cross section in a naive 'fully parton-level' way omitting the proton 
at the top of the diagram because then one would have three particles in the initial state 
(whereas the standard framework for calculating a cross section requires two particles in 
the initial state). Furthermore, by deleting the proton at the top of the diagram one would 
then be neglecting the important fact that the two partons on this side are tied together in 
the same proton (as was pointed out in [19]). The use of proton wavefunctions or hadronic 
amplitudes in the calculation of DPS-type graphs was discussed long ago in [27], and has 
been discussed more recently in [6, 19]. We utilise the approach and notation of [27] in our 
work. That is, we assign a wavefunction factor (p to the p — t- qqX vertex that is assumed 
to be strongly damped for values of the parton transverse momentum and virtuality larger 
than the hadronic scale A ~ 1/Rp ■ In our case the factor (p is a matrix in spinor space, 
and also carries a label x that describes the spins of all of the particles in X. 

In the following we will take a number of steps to simplify the calculation as much 
as possible. First, we will largely ignore considerations of colour, and will suppress colour 
indices, factors and sums where they appear, in order to avoid the proliferation of too many 
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indices. Second, we will take take the four-momenta squared of the two off-shell photons 
to be the same, and refer to this common four- momentum squared as . Finally, we will 
take the colliding protons to be unpolarised, as is the case for the colliding protons at the 
LHC. 

As in [20], we apply a Sudakov decomposition to all four- vectors used - that is, we 
write an arbitrary four- vector Y in terms of two lightlike vectors p = -^(1,0,0,1) and 
n = -^(1, 0, 0, —1) and a transverse part denoted as V: 

y = y+p + T/-n + V (2.1) 

Rather than proceeding to calculate the cross section contribution from figure 4(a) di- 
rectly, we instead begin by calculating the cross section contribution a2v2 associated with 
the Feynman diagram in figure 4(b). In this diagram, two nonperturbatively generated 
quark-antiquark pairs produced by colliding protons interact via two separate qq — )• 7* 
hard processes. We should be able to express the cross section of this process in terms of 
'nonperturbatively generated parton pair' r-space 2pGPDs r(xi,X2; A) and hard subpro- 
cess cross sections a as follows: 

o"2^2(s)=y" dxxdx2dyxdy2aqq^^'.{s = XYl)\S)aqq^^*{s = X2y2S) (2.2) 

-j^—:^Tp{xx,X2\^)Vp[yx,y2\ -A) 

Helicity labels are omitted in the above schematic expression, but they will be included 
in the full calculation below. By using the fact that the expression for the cross section 
must end up in this form, we can establish the connection between the vertex factor Lp and 
the 'nonperturbatively generated parton pair' r-space 2pGPD F. We shall need to make 
use of this relationship when we come to study figure 4(a). 

Note that a calculation of a2v2 has already been performed by Paver and Treleani in 
[27] for the case of spinless partons, and by Mekhfi [28] and Diehl, Ostermeier and Schafer 
[5, 6] for the case of partons with spin. We follow closely the approach of Paver and 
Treleani, and our calculation of o"2„2 can be considered as a brief review of the method in 
[27]. 

We will neglect the proton mass with respect to the total centre of mass energy 
and work in a frame in which A is proportional to p, whilst B is proportional to n, A = 
A+p, B = B^n. One can directly write down the following expression for the cross section 
contribution from figure 4(b), (72v2{s)'- 

a2v2is) = 2(2l)io Yl I d'^Ad^Bd^Jid'^J25^^\A + B + J1 + J2-A- B)5{JI - Q^) 

X <5( J| - g2);\^X7Mi/^2(^^ A, B, Ji, J2)M''"'^'^'' {A, B; A, B, Ji, J2)* (2.3) 

where: 

M^"'>''^'^A,B;A,B,Ji,J2) (2.4) 

Tr [r/^i(Ji)^,ypg(ai,a2,i)^,r/^^(J2)fe(/9g(fe2,bi,^)l^i] 
(27r)4 D{ai)Dia2)D{h)Dib2) 
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D(a) =a + ie, 



02 



A — A — ai bi = Ji — ai h2 = B — B + ai — Ji 



(2.5) 



(2.6) 



The vertex factors 99 ensure that the quark and antiquark hues with momenta Oj and 
hi have small virtuality. Given that this is the case, we can rewrite the slashed vectors in 
(2.4) as sums over outer products of particle or antiparticle spinors (as appropriate), using 
the completeness relations. Then we have: 



M^^^^^'^A^B-A.B.Ji,,]- 
-74 



(2.7) 



(02^2 



J2 



u''^{ai)ip^{ai,a2,A)v''^{a2) 



u'^{b2)ipUb2,bi,By^{h] 
D{bi)D{b2) 



D{ai)D{a2) 

The Si and ti are quark or antiquark helicity labels, and the Mqq^^* factors are 'hard' 
qq — )• 7* matrix elements. The hard matrix elements should be evaluated with initial state 
partons having small (i.e. hadron scale) transverse momenta and off-shellness - however, we 
make the approximation in the matrix elements that the initial-state partons are on-shell 
and collinear, which only corresponds to a small relative error 0{h?' /Q'^) <^ 1. 

Consider now the integrations over the longitudinal parts of ai - i.e. and . It 
is not hard to show that the integration over is restricted to values of order A^/Q by 
the vertex factor {p{ai,a2, A), whilst (p{bi,b2, B), D{bi), and the Mqq-^^* are practically 
constant in this range (and approximately equal to their values with set to zero). 
Similarly, the integration over is restricted to values differing from by ~ A^/Q by 
the vertex factor </?(6i, 62, B), with (p{ai, 02, ^4), D{ai) and the Mqq^^* being approximately 
constant and equal to their values at = in this range. This allows us to write: 



M 



(A, B- A, B, Ji, J2) ^ Kf^r ('^ip^ -^r^ ^ -^i) 

cPa\ \ f da^ u^^{ai)ipp{ai,a2,A)v^^{a2) 



J2 



(2vr)2 

dbt u'^{b2)v;{b2M,By{biy 



27r 



D{bi)D{h2 



2tt 



--Ji 



D{ai)D{a2) 



(2.8) 



Define: 



da^ n^i(ai)(^^(ai,a2, A)f*2(a2) 
D{ai)D{a2) 

Then we can write M'^'^^'^^'^[A, B; A, B, Ji, J2) in a more compact form: 
M^^^''''^A,B;A,B,Ji,J2) 

^T.Ki^r^JtP^Jin ^ Ji)M'^i^';i{4p,J2n ^ J2) 

d'^a-i 



(2. 



(2.10) 



(2vr) 



2 V'^-g'f ^2^, ai, «2, i")Vg*f (J2-, Jf, &2, bi, ^+) 
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We now insert (2.10) into (2.3), and make use of the following relation in the resulting 
expression: 

Mlf^t;i{Jtp, J{n ^ J,)M;li^>^'{Jtp, J{n ^ Ji)(2vr)<5( - Q^) (2.11) 

= a^^!^^'*^^^^(s = 2J+Jr)4J^Jr 

^qq'^^^*'^^'^^ is the qq — )• 7* 'cross section' with q, q, 7* helicities si, ii, /ii in the matrix 
element, and in the conjugate matrix element (note that if si 7^ si and/or ti 7^ ii 

this is not a cross section in the strict sense). (5"^^!!^^^*^'*^''^^ is related to the spin-averaged 
qq — )• 7* cross section Uqq-^^* by: 



- si,ti;si,ti;/^l 



a, 



qq—i-'y* 



2(Tgg^^*5sj -t^J^^ (2.12) 



The result of inserting (2.10) into (2.3) is: 
^2.2(^) = ^^^^ E j d'Ad'Bd'JiKy^f''''''{s = 2J^ (2.13) 

SiUsiliXy 

^ '^qq^'Y* \^ — ^'^2 ■-'2 )^'^2 "^2 J (271)"^ '^''^^ i"l)"'2)^ ) 

><^ql'^i''''''^'is = 2J+J^) j d^aid^bid^Ad^AdA-d^BdB+ 
E ('^1'' 4^ ai, a2, i-)^;Sf ^(J+, 4,a^ + A, - A, 1-)4J+ 

X 

E ^pm'{J2^JiMM,B+)^pf^f{j:,,4M + A, bi - A, ^+)4Jf J^^i?" 



X 

X 

X 

7 



In the second line of (2.13), we have introduced the transverse variable A via ai = 
ai + A, and converted the integral over J\ to an integral over h\ using ai + bi = Ji. 
We have also made use of the fact that s = 2A'^B~ . Let us define the 'nonperturbatively 
generated parton pair' r-space 2pGPD according to: 



ei?''^^^ f ^A:A^ / dA~d^Ad^mrp^^^^{J+,J+,a^,a2,A-,A) 



X 



(2.14) 



X i^l-^t^'i^t, 4, «i + A, a2 - A, , A)J+4A+ 



We also introduce the following scaling variables: 

xi = J+/A+ X2 = 4/A+ yi = 4/B~ y2 = 4/B- (2.15) 
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Changing variables in (2.13) to the scahng variables, replacing appropriate combina- 
tions of ^ps by Fs according to (2.14), and using the obvious relation TpjH'*^*^ (xi, X2; A) = 
rpq^'^'^^^{x2,xi;—A.), we finally obtain: 

cr2v2{s) = ^ j dxidx2dyidy2aq\'^:ll^'^^'''\s = xiyis)aql'^^^^ (2.16) 

X I ^r;|,1'^'^^"^(xi,X2;A)r;^|f *^(yi,y2;-A) 

The cross section is of the anticipated form (2.2). The most important result of this 
preliminary calculation is the definition of the 'nonperturbatively generated parton pair' 
r-space 2pGPD (2.14), which we shall make use of later. 

The calculation of the cross section contribution associated with figure 4(a), cr2vi{s), 
proceeds in a very similar manner to the calculation of a2v2is)- Once again we work in 
a frame in which A = A^p and B = B'^n. We can directly write down the following 
expression for the cross section: 

c72vi{s) =^(^ E / d''Ad'j,d'j2S{jf - Q^)5{4 - Q2)5(4)(i + j^ + j^_A-B) 

(2.17) 

X A^^;xw/^2 {A, B; A, Ji, J2)M^'^'''^'^ (A, B; A, Ji , J2)* 

where: 

7W^'>^^i^2(^,S;i, Ji, J2) (2.18) 
^ / ^^'^K^i9'^(ai,«2,^)^2^^''^'''n«2aii?^ JiJ2))/[I)(ai)Z)(a2)], 



r^'^i^2(a2aiS ^ J1J2) = iHeQq^g: 



f^^{j2)hh(B)kr,M^ 

Dih)D{b2) 



(2.19) 



The lines with momentum ai are restricted to small virtuality by cpA, so we can de- 
compose the slashed a, vectors in (2.18) into outer products of particle or antiparticle 
spinors: 



M^'^''^f^{A,B;A,Ji,J2) 

d'^ai \ 12*1(01)99^(01,02,^)^^2(02) 



D{ai)D{a2) 



(2.20) 



^s2SiA;/ii^2 ^Q^oii? — )• J1J2) is the matrix element for qqg — )■ 7*7* with initial quark and 
antiquark having small transverse momentum and virtuality. 

For reasons similar to those leading to equation (2.8), we can move the integration 
such that it only acts on the part of (2.20) in square brackets, and set a||~ = in the rest of 
the integrand. Provided that that ^ A^, we can perform an analogous operation for 
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the tti integration. The reason for this is that when Ji^ ^ A^, the transverse momenta of 
the Oj hnes (constrained to be of order A by if a) are neghgible compared to the transverse 
momenta of the bi and J, hnes in J^ '^2SiX;fiifi2 (^^^aiB — t- J1J2), so we make only a smah error 
by setting to zero in this factor provided Ji^ 3> A^. Applying these approximations: 



M^'^'''''^A,B;A,Ji,J2) 

Si •' 



(2.21) 



(27r)2 2tt 



D{ai)D{a2) 



We identify the final factor in (2.21) as the integral of ^lJp over ai. Writing out the 
denominator factors in Ad^'^^^^'^^'^'^ {a2aiB — )• J1J2) explicitly we have: 



M^'^''^^'^{A,B;A,Ji,J2) 



(2.22) 



da^ 
'2^ 



d'a^/{2TTYrp"''''{at,o+, ar,a2, A^ 
[2(J+ - a+) Jf - + ie][2{at - J+)J^ - Ji^ + ie] 



where: 



r^2srX-M^^^^^B ^ J1J2) ^ i5^,(eQ,)V2(a2)^;^(J2)l^2^A(^)^i^;i('^i)^'Hai) (2.23) 

Examination of the denominator factors in (2.21) reveals that the majority of the 
contribution to the af integration comes from the region af ~ . For this reason we can 
set af = in the numerator before evaluating the af integral using contour integration: 

M^''^'''^'^{A,B;A,Ji,J2) (2.24) 

^^ilf d2ai/(2vr)Vp^^^^(J+,J2+,ai,a2,^-) 
s^ Li 

r-^-^^^^^^^(a2aii?^ JiJ2)L-=o,,,=o,,+=^+ 
2(Jf + J2")Ji2 

We are interested in the behaviour of M^''^''^^''^ {A, B; A, Ji, J2) when Ji^ <C (but 
still A^) such that all of the internal particles have transverse momenta and virtualities 
much less than Q. In this limit we can use spinor completeness relations to split T up 
into two qq — )• 7* matrix elements and one g ^ qq matrix element, with the quark and 
antiquark having small transverse momenta and virtuality 0(1 Ji|) in each matrix element. 
The quark and antiquark transverse momenta and virtualities can be set to zero in the 
'hard' qq — )• 7* matrix elements with only a small accompanying error 0{J\^ /Q"^), but we 
must keep the term proportional to J\ in the g ^ qq matrix element as this vanishes in 
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the limit J\ — )• 0: 

M^''^^'^'^{A,B;A,Ji,J2) (2.25) 
-i [/(i2ai/(27r)Vp'''''(a^4,ai,a2,i-)] M^-^%'^{B- J{n + Ji,J^n + J2) 



Sit 2 



2(Jf + Ji^ 



Having inserted (2.25) into (2.17), we use (2.11) and the following connection between 
Mg^qq and helicity-dependent unregularised splitting functions in the result [29]: 

-A^^4f (S; Jfn + Ji, J^n + J2)Mf^l''' (B; J{n + Ji, J^n + J2) (2.26) 



9„2 pX—>-t2tl ,t2ti 



Jo 



This yields: 



0-21,1 (s) = ^ 



X 
X 



qq—^y* 



(27r)i2s 



S = 2J+Jf )J+(T, 



2 pA^t2tl,t2tl 



Jo 



S — J2 )J2 

(2.27) 

1 



(2^)3^+2^- 



si,ii;Si,ti;/^l 



qq—^y* 



(s — 2J2 J2 



(2.28) 



^^'^^(^1+, J2^,ai + A,a2 - A,^-)J+J2+^+ 



2 pX^htiMi 
iJs ^g^qq 



Jo 



Jl + J2 / J^ 



In (2.28) we have once again introduced the transverse variable A via the same relation 
as in the 2v2 case. We recognise the object in square brackets in (2.28) as the integral of 
the nonperturbatively generated parton pair r-space 2pGPD over A. If we make a change 
of longitudinal integration variables in (2.28) to the scaling variables (2.15), then we finally 
obtain: 

0-21,1 (s) = ^ dxidx2dyidy2afq'!^i^l^'^^''^^ {s = xiyis)aql'^^^^^ (2.29) 



^i^i ti 



(2vr) 



_pSl 



2 pm 



{xi,X2; A) 



^Pg^qq (?/2) S{1 - y, - ^2) J — 



We have restricted our integration over Ji to the range 
corresponds to the range over which our approximate expression for the matrix element 
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(2.25) is valid. The contributions to a2vi coming from Ji^ values outside this range do not 
have the same 1/Ji^ structure. 

The integral over Ji in (2.29) gives rise to a large transverse momentum logarithm 
log(Q^/A2), whilst the integral over A gives a prefactor of order ~ (since the 

nonperturbatively generated parton pair r-space 2pGPD only has support for momentum 
imbalances of order Aqcd)- Thus, as we asserted at the beginning of this section, there is 
a part of the cross section expression for figure 4(a) that is proportional to log{Q'^ / A"^) / Rp 
and should be included in the LO DPS cross section. 

Note that the quantity / d^AFp^*^'*^*^ (xi,X2; A) /(27r)2 is equal to the b-space non- 
perturbatively generated parton pair 2pGPD evaluated at zero transverse separation, 
-psiS2,siS2 b = 0). This appears to indicate that the 2vl contribution to DPS probes 

nonperturbatively generated parton pair 2pGPDs at zero parton separation. In fact, the 
result (2.29) actually corresponds to a broad logarithmic integral over values of that 
are < but > The 6-space 2pGPD evaluated at 6 = appears in (2.29) because 

the r-space nonperturbatively generated parton pair 2pGPD dies off rapidly for A^ ^ A^, 
which is equivalent to the 6-space nonperturbatively generated parton pair 2pGPD not 
containing any fluctuations with length scales <C Rp. Then we can approximate Tp{b) for 
the relevant values of b in (2.29) by Tp{b = 0). 

If one assumes that diagrams of the form of figure 1(b) are the only diagrams of the 
'2vl' type that contribute to the DPS cross section at leading logarithmic order, then a 
generalisation of the result in (2.29) yields the expression below for the contribution of 2vl 
graphs to the LO DPS cross section^ : 

=2 X T E ^ dk J dxidx,dy,dy,-r-r-r^ (2.30) 

X ^iiji^Ais = xiyis)ai2j2^B{s = X2y2s) 

X ^H^^P,^,,,, (^) ( ^; k\Q^) ( Sl; k^Q'- 

y[ + y'2 ■'^■'^ \yi + y2J ^^Kvi J \y2 

(^;A^Q^)i^- (g;A^Q^)^;;:L.,(x;,4,^>=o;A^) 

(x; k'^,Q'^^ are the Green's functions of the DGLAP equations - i.e. a set of func- 
tions obeying the DGLAP equations with the initial condition Dj (x; k"^, k'^) = 6ij 5(1 — x). 
Pp *^^gp(x'^, X2; 6 = 0,A2) represents a nonperturbative initial condition for the two inde- 
pendent ladders in figure 1(b). In (2.30) we have re-inserted the symmetry factor m/2 that 
has been omitted in earlier discussion in this section (m = 1 if the two hard processes are 



^Note that here and in the rest of this section we will take the scales associated with the two hard 
scales to be equal, Q\ = Q% = Q^. We will comment in section 4 on the generalisation of the results 
of this section to the case of unequal scales. Note also that we only write down the unpolarised diagonal 
contribution in colour, flavour and spin space here. The contributions associated with spin polarisation 
(either longitudinal or transverse) and flavour interference are expected to have a similar structure. On the 
other hand, it is known that the colour correlation/interference and parton type interference contributions 
will be suppressed by Sudakov factors, as is discussed in [5, 6, 26, 30]. 
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identical, and m = 2 otherwise). There is an additional prefactor of 2 in (2.30) because 
there are two sets of 2vl graphs that give equivalent contributions - in one set the non- 
perturbatively generated parton pair emerges from the 'left' proton, whilst in the other it 
emerges from the 'right' proton. 

Equation (2.30) can be written in a more compact fashion as: 

(^fA^){s) =2 X y j dxidx2dyidy2CFi^j^^A{s = xiyis)ai^j^^B{s = X2y2s) (2.31) 

xDi^'Hyi,y2;Q') J ^rif^,epi^,,x2,A;Q') 

=2 X y ^ j dxidx2dyidy2ai^ji^A{s = xiyis)ai2j2^B{s = X2y2s) 
X bi^^^{yuy2;Q^)V';^^,,p{x^,x2,b = Q-Q^) 



where: 



'<r;:;Uw,4,6;A^) 

=E r ..^^##ii«±^ (2.33, 

As mentioned in section 1, and as will be explored in detail in section 3, there are 
additional diagrams of the '2vl' type that contribute at leading logarithmic order to the 
DPS cross section, aside from those represented by figure 1(b). These involve crosstalk 
interactions between the two nonperturbatively generated ladders. Equation (2.30) (or 
(2.31)) therefore represents only part of the 2vl contribution to the LO DPS cross section. 
For the moment, however, we'll limit our discussion to just this part. 

A necessary requirement for (2.30) (or (2.31)) to be valid (at least as an incomplete 
part of a contribution to the DPS cross section) is that the independent two-ladder 2pGPD 
Pp *^^gp(xi, X2; b, Q"^) should be smooth on distance scales ^ i?p ~ 1/A (or equivalently that 
the corresponding distribution in terms of the transverse momentum imbalance A is cut off 
at values of order A). This appears to be a somewhat reasonable requirement - at the scale 
A there is only this scale available to set the size of the A profile for ^^^^^(xi, X2; A, A^), 
and the evolution equation for independent two-ladder 2pGPD (which is just the 'double 
DGLAP' equation of [21, 22] with the 'single PDF feed term' defined in [31] removed) 
preserves the transverse profile. In any case, such behaviour for P^ ^^^^^(xi, 2:2; A, Q^) 
would appear to be required in order to get the necessary prefactor of order 1/-Rp in the 
2v2 contribution to DPS, which is calculated according to the following expression (for the 
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diagonal unpolarised contribution): 

'^{A%^)(^) =Y j dxidx2dyidy2cri^j^^A{s = xiyis)ai^j^^B{s = X2y2s) (2.34) 

X j (^r;^^^dep(^l,^2,A;Q2)r^l^^^^^(y,^ 

Y ^ y dxidx2dyidy2cri^j^-,A{s = xiyis)ai^j^^B{s = X2y2s) 
X / d'bTl%,^^{x,,X2,b;Q')Tl%,^^^^^^ 

If one assumes that Tp j„rfep(^i' ^2; &i Q^) can be factorised into a longitudinal piece 
l)p j^^gp(xi, X2; Q^) and a flavour-independent transverse piece F(&), where F{b) is a smooth 
function of radius Rp normalised to 1, then (2.30) and (2.34) become: 

^f^fm(^) ^- / dxidx2dyidy2d-i^j^-^A{s = xiyis) (2.35) 

V , 1 I ^ a-eff,2v2 J 



X 



^i2j2^B{s = X2y2s)Dl'}laep(^uX2; Q^)Dl]indepiyi^ V^'^ 



D,2vl 



(s) =2 X ^ / dxidx2dyidy2ai,j,^A{s = xiyis) (2.36) 

^ 7""^ '^eff, 2vl J 

X ai^j,^B{s = x2y2s)Dp^Hyi,y2;Q^)Dl%dep(^i,x2;Q^) 



where: 



^ EE / d'b[F{b)]' = I -^[F(A)]2 (2.37) 
/,2«2 J J (2vr)^ 

^F{b = 0)= [ ^[FiA)] (2.38) 



-F(A) is the Fourier transform of F{b). We see that the geometrical prefactors for 
the two different contributions to the DPS cross section are different in general, (Jeff, 2v2 7^ 
(7eff,2vi- If one assumes that two nonperturbatively generated ladders are to some degree 
uncorrelated in transverse space, F{b) is given by a convolution of an azimuthally symmet- 
ric transverse parton density in the proton p{r) with itself, where p{r) must be normalised 
to 1 in order to ensure the appropriate normalisation of F{b): 

F{b) = j d^rp{r)p{b - r) (2.39) 

Then, if one takes the Gaussian form exp[— r^/(2i?^)]/(7ri?^) for p (with R a constant 
parameter), one finds that cre//^2vi = o'e//,2D2/2 - that is, the 2vl contribution receives a 
factor of 2 enhancement over the 2v2 contribution from the geometrical prefactor alone 
(in the next section, we'll discover that the 2vl contribution is further enhanced at low 
X as a result of the crosstalk interactions on the two-ladder side that are allowed for this 
contribution). The ratio CTef f,2v2/ <^ef f,2vi does not depend much on the precise shape of p 
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- for example, one obtains 2.18 if p is a top hat :^^Q{R — r), 2.32 if p is the projection of 
an exponential J dz g^j^a exp(Vr2 + 22/i?), and 1.94 if p is the projection of a hard sphere 

— r2/i?2)^/2 0(i? — r) (with R once again a constant parameter in these expressions). 
It is important to bear in mind, however, that in order to obtain an enhancement that is 
roughly a factor of 2 we have had to make a number of assumptions whose validity is 
somewhat uncertain (this is particularly the case for the assumption (2.39)). There could 
be some 'clustering' of the nonperturbative partons in transverse space, which would tend 
to increase (yeff,2v2/ (^eff, 2vi- Alternatively it is not inconceivable that the probability to 
find two nonperturbative partons separated by small distances <^ Rp could be smaller 
than the probability to find them separated by distances of order Rp - in this scenario 

/o"e//,2i)i would be reduced. 

3 Crosstalk between Ladders in the 2vl Contribution 

In the previous section we demonstrated that there is a leading logarithmic contribution 
to the DPS cross section associated with diagrams in which a single parton ladder from 
one proton splits into two, and then the two daughter ladders interact with two indepen- 
dent ladders from the other proton (that are only connected to one another via low-scale 
nonperturbative interactions). It is suggested in a number of works [19, 25, 32] that these 
diagrams are the only ones involving a single 1 — )• 2 ladder branching that give rise to a 
leading logarithmic contribution to DPS. Here, we show that there is also a leading log- 
arithmic contribution to the DPS cross section associated with diagrams such as those 
in figure 3 in which the two nonperturbatively generated ladders talk to one another by 
exchanging partons, provided that the crosstalk occurs at a lower scale than the scale of 
the 1 — ?■ 2 ladder branching. There are two types of crosstalk that are possible, which 
are illustrated in the simple diagrams in figure 5(a) and (b) - we'll call these off-diagonal 
real emission and virtual exchange processes respectively. As in the previous section, we'll 
demonstrate that there is a leading logarithmic contribution from diagrams such as figure 3 
by examining one of the simplest possible diagrams of the appropriate type - namely, that 
of figure 5(a). We will find that there is a large DGLAP logarithm associated with both 
the 1 — )• 2 splitting and the off-diagonal real emission ('crosstalk') processes in the figure, 
and that this is associated exclusively with the region of integration in which the partonic 
products of the off-diagonal real emission have much smaller transverse momentum than 
the products of the 1 — )• 2 splitting (and all of these transverse momenta are 3> but 
« Q2). 

In our calculation, we'll ignore considerations of colour for simplicity, just as we did 
in section 2. However, the colour structure of crosstalk processes is quite nontrivial, and 
is important when considering the size of such contributions to cross sections. The colour 
structure of crosstalk processes has been considered previously in the context of twist-4 
contributions to DIS in [33-36], and in the context of DPS in [24, 26, 37]. We will make 
some comments with regards to the colour structure of the crosstalk processes at the end 
of this section. 
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Figure 5. (a) Simple 2vl diagram including an 'off-diagonal real emission' process, (b) Simple 2vl 
diagram including a 'virtual exchange' process. 



As in section 2, we work in a frame in which A = A'^p, B = B n. The cross section 
expression associated with figure 5(a) is: 

cTXTis) = ^^^lyo^ j d^Ad^rd^Jid^J26{j! - Q^)6{Jl - Q^)6{r^) (3.1) 
Si^)^A + B-A-r- Ji- J2)Ml''^'^'^'''iA, B;Ji,J2,r,A) 

where: 

M^^^'^'^^'^{A,B;Ji,J2,r,A) (3.2) 



(27r)4 D{ai)D{bi)D{a'2)D{b2)D{a2) 



A^i^''^™^(Ai?;Ji,J2,r,A) (3.3) 



,^(B)^i^;^(Ji)^i^;^(r)^;(^^(5i,a2,A)^24(J2)^2' 



(27r)4 D{~ai)D{bi)D{d[)D{b2)D{~a2) 



Following a procedure that is similar to that leading to equation (2.22), and is valid in 
the region of transverse momentum integration in which jf, J^i''^ ^ (c"^ equivalently 
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ttg, ^2 ^ ^ can write down the following approximate expression for A4l- 



(3.4) 



Sis, 



'2^ 



d'ai/{2TTyiP;''^''{at ,a+ ,ai, a'^, A' 



■7-^2 

'L 



2SiA;/xl^2M3 / / D , T T M 



[2(J+ - at)Ji -Ji + ie][2{a+ - J{){B- - Jf ) - Jf + ie] 
1 



[2(J+ + J2+ - a|)(-r-) - r2 + ie] 

where TL{a2aiB — JiJ2r) includes all of the numerator structure of the ML{a2aiB — )• 
JiJ2f') matrix element. 

Performing the integral using contour methods, and making use of the fact that 
the overall integrand is strongly peaked near = whilst the numerator factor 71 is a 
relatively smooth function in this region, we obtain: 



(3.5) 



-ir 



S1S2 



dPa^/{2TTf^;''^''{Jt. 4+, ai, 4, i-) 



_S2SlA;/ll/X2/13 / / D , T T M 



In the region of integration in which j\ <^ Q^, we can drop the second term in the 
denominator factor [J^ + -^pr + r+]. Also, when jf, J^^r'^ <C Q"^, we can approximately 
decompose 71 as follows: 



(3.6) 



Mf^'^i;^[4p. J2n ^ J2)M^''' {a'+p,4p - r, (, 



2-5>S2/i3 / ^ + , 



Oq"^ — at 



)p + r) 



X A^J4f {B; J^n + J^,J-n-Ji] 



Performing a similar sequence of operations for Mji, we obtain the expression below 
that is vahd for A^ < af , b^, 63 < Q^, or equivalently A^ < jf , j|, <C Q^. 



A^^'^™^(Ai3;-/i,J2,r,^) 



d2a;/(27r)^V;'^"''(«i+, J2^ a'l, «2, ^" 



(3.7) 



A^fir.r^ (a'l+P, a+p - r, (a'+ - a+)p + r) 



2J," 



X >[™(J+p, Jfn ^ Ji)>[;™(J2+p, J^-n ^ J2) 
X (S; Jfn - J2, J2~n + J2) 
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Given that the transverse momenta of the partons emerging from the g ^ qq branching 
process are different on the left and right hand sides of the cut in figure 5(a) (=b Ji and 
=F J2 respectively), we will require a generalised version of the relation (2.26), which reads: 

'^^ > -Mg49*' -^1^ + ^1' -^2^ - Ji)M;^gl''' {B- J^n - J2, J^n + J2) (3.8) 



Note that in the off-diagonal emission process, the partons emitting the gluon in the 
amplitude and conjugate do not in general have the same plus momentum (and indeed 
are not of the same type). This means that the product of Mq-^qg and Mq^qg from 
the left and right hand sides of the diagram does not give rise to a conventional splitting 
function multiplied by the appropriate transverse momentum squared, as occurred in (2.26). 
Instead, one obtains: 




'-^Ml^^';'^ {a'+p, atp - r, (a'+ - a+)p + r) (3.9) 
xMll^;'''' {d'+p, afp - r, {~a'+ - at)p + r) 

9«2T/sis2-^siS2;At3 / 1 ».2 

^9s ^I,q^q \A+' A+' A+ I 



°'i4-!>si.52;M3 / a+ 



where Vj]^^^ ' y^V^ represents some kind of generalised splitting function, 

that satisfies the following relation: 

Furthermore, since the partons emerging from the hadronic blob in figure 5(a) do not 
in general carry the same momentum on the left and right hand sides of the diagram, the 
process in figure 5(a) probes a two-parton PDF that is not diagonal in x. It is defined 
according to: 

- (^ E / dA-d^Ad^a,d^a',^ata',^d[^dtA+ (3.11) 

X i^p-q^^iaf, a'+,ai,a'2, A- , A)^P*p''^^'^^ {d[+ , d^ , a[,d2,A-,A) 

Note that this distribution is somewhat similar to the four-quark matrix element that 
is probed in the twist-four contribution to Drell-Yan, and that is defined in [38-41]. Here, 
however, we do not absorb two powers of the strong coupling constant gg into the four quark 
matrix element, as is done (and makes sense) in the context of the twist-four contribution 
to Drell-Yan. 

Inserting (3.5), (3.6), and (3.7) into (3.1), and making use of (3.8), (3.9) and (3.11), 
we find that the contribution to axT coming from the region of transverse momentum 
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integration with <C r^, jf , <C is: 

crxris) = ^ j dxidx2dyidy2&qg^:^l^'*^'^^{s = xiyis)ag^q'^^^^ 

(3.12) 

^ J ^^'i^/,g4g ^ ^'^^(^1) ^'i) ^2)rp|g|' ^ ^ (3^1, 2:2, x'l) 

(2/2) 6(1 - 2/1 - y2)J y dJi dr ^2j2(^^+^)2 

In the region of transverse momentum integration in which ^ 17^,^2) the transverse 
momentum integrand simplifies as below, and we obtain two large DGLAP logarithms from 
this region: 

2 Jle^(Jl+r) fQ' dJ^ f-"^" dr-" _^ ^ f 



dJ^'dr' \ X'r o > / ^ / ^ = log' ^ (3.13) 

Two large DGLAP logarithms implies a leading logarithmic contribution, since there 
are two powers of in (3.12). Thus, there is a leading logarithmic contribution to the 
DPS cross section coming from the region of figure 5(a) in which r"^ <^ j\ (i.e. in which the 
scale of the off-diagonal real emission process is strictly smaller than the scale of the 1 — )• 2 
branching process). Note that it is only this region of transverse momentum integration 
that gives rise to a leading double logarithm - other regions only give rise to either a single 
logarithm, or no logarithm at all. The single DGLAP logarithm is essentially associated 
with a logarithmic integral over r only, and this should be absorbed into the four-quark 
matrix element in the 'conventional' twist-4 contribution to double Drell-Yan. 

Aside from the process in figure 5(a) involving an off-diagonal real emission, the process 
in figure 5(b) involving a virtual exchange also gives rise to a leading double logarithm, 
provided once again that the virtual exchange process occurs at a lower scale than the 
1—7-2 branching. This is straightforward to show using a procedure similar to the one 
we have used above. Generalising these results, we find that in the most general 2vl DPS 
diagram, all possible types of parton exchange are allowed inside the two ladders emerging 
from one of the protons at leading logarithmic order, provided that they occur at a lower 
scale than the 1—7-2 ladder branching occurring in the other proton. Schematically, the 
LO (diagonal unpolarised) cross section expression for the 2vl contribution to DPS is: 

D,2vi, . n^'"'^ \^ Z"*^' wi.2"s(fc^) dx[ dx'2 dy[ dy'2 

V,B)(«)=2x^ ^ J^^ dk ^ J dx,dx2dy,dy2 — — — — (3.14) 

liijii'ij'i 

X <5"nji->A(s = xiyis)ai2j2~^B{s = X2y2s) 

2/1+2/2 V 2/1 +2/2/ ^1 V2/1 / ^2 \^2/2 

X ^1 kx'i,4;^'i,A.2) 
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T^i^^^{xi,X2',xi, fj,'^) is a four-parton matrix element whose evolution involves all possible 
exchanges between these partons in an axial gauge^ [42]. Dj (x;/c^,Q^) are the Green's 
functions of the DGLAP equations - i.e. a set of functions obeying the DGLAP equations 
with the initial condition Dj (x; /c^, /c^) = 6ij 5{1 — x). 

Note that at the next-to-leading logarithmic (or NLO) level, one needs to append an 
extra term to (3.14) that is of the following form: 

j dxidx2dxidyD'l{y, Q^)T'i{xi,X2;xi,Q'^)aijk^AB{xi,X2,xi,y) (3.15) 

This is essentially the 'conventional' twist-4 contribution to the pp — t- AB + X pro- 
duction cross section. At the level of total cross sections, the DPS contribution to the 
production of AB cannot be distinguished from the conventional twist-4 contribution, and 
the two should really just be considered together as components of the 0{h? /Q"^) correction 
to the pp — >• AB + X cross section. 

Let us now discuss the issue of colour in the evolution of the four-parton (twist-4) 

i' i' 9 

matrix element F^^ {xi,X2'-,xi, ). We recall that, for the 2pGPD with finite fo, every 
distribution that does not have the partons with the same light-cone momentum fractions 
on either side of the cut paired up into colour singlets is suppressed by a Sudakov factor - see 
[5, 6, 26, 30]. This factor arises in axial gauge because there is an incomplete cancellation 
of the soft gluon region between (diagonal) real emission diagrams and virtual self-energy 
corrections in the colour interference/correlation distributions [30]. In physical terms, it 
occurs because such distributions involve a movement of colour by the large transverse 
distance h in the hadron [26]. 

i' i' 

In the twist-4 matrix element T^^ ^{xi,X2]Xi, ^ ) there is no such Sudakov suppres- 
sion of colour interference/correlation distributions. The extra diagrams that are allowed 
in the evolution of this distribution (i.e. the off-diagonal emission and virtual exchange 
diagrams in axial gauge) provide extra soft-gluon divergences that cancel any remaining 
divergence from adding the diagonal real emission and virtual self-energy diagrams to- 
gether. The soft divergence in both real emission diagrams (diagonal and off-diagonal) is 
positive, whilst that in both virtual diagrams (self-energy and exchange) is negative, and 
in the sum the positive and negative contributions always cancel each other out. We can 
see why this cancellation occurs physically as follows. In the operator definition of the 
twist-4 matrix element, the four operators corresponding to the partons all lie on the same 
lightlike line, with no transverse separation between any of them. Note that this does not 
correspond exactly to the physical situation - in (3.14) the transverse separation of the 
partons in F^^*^ (x'j^, x'j^, k"^) is in fact more like l/|fc| - but for the purposes of obtaining 
F^^*^(x']^, x'^, A;^) by solving the evolution equation at scales fj? < k"^, the distinction is 

i' i' 

unimportant. The fact that the four operators/partons in F^^ ^(xi, X2; xi, ) are on top of 
(or at least very close to) one another in transverse space means that soft longwave gluons 

^In a covariant gauge, such as Feynman gauge, there are further diagrams that contribute to the evolution 
due to the presence of a nontrivial Wilson line in the definition of the operator. These diagrams involve 
gluon connections to the Wilson line. 
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Figure 6. A process that can bring about a colour recombination in the four gluon state. On the 
diagram we have indicated the colour flow in the large Nc limit. 

can only resolve the total colour of all of them. But the summed colour of the four partons 
must be zero, since the proton is a colour singlet object - therefore the effects of soft gluons 
must cancel, as is indeed observed in practical calculations. The cancellation of soft gluon 
divergences in the twist-4 matrix elements has been discussed before, in [33, 43, 44] (for 
example) . 

Let us now focus our attention on the region of small x (which is perhaps the most 
relevant region of x for DPS processes at the LHC). It is well known that in this region 
the gluons dominate, so we will only consider these partons in what follows. We have seen 
that the colour correlated/interference twist-4 distributions are not Sudakov suppressed - 
however, in the low x region, the distributions in which two pairs of gluons are in colour 
singlet configurations tend to win out. This is because the colour factors in the anomalous 
dimensions for these distributions are larger (see section 3.2 of [33] or section 5.1.3 of 
[6]). Bear in mind that in figure 3 at scale k'^ the nonperturbatively generated partons 
with identical x fractions must be in a colour singlet state if one wants to avoid any 
Sudakov suppression. By combining two off-diagonal real emission processes put together 
with two diagonal real emission processes as in figure 6, it is possible to achieve a 'colour 
recombination' on the two ladder side at scales lower than fc^, and alter the way in which 
the parton legs are grouped into two sets of colour singlets. For example, in figure 6 the 
grouping is changed from (14)(23) afterwards to (13)(24) before, using the leg labelling 
conventions from the figure. Such a colour recombination is not disfavoured from the point 
of view of evolution before and after the process - however, it is itself suppressed by a 
colour factor equal to 1/{Nq — 1) [33-37]. This colour factor suppression is associated with 
the fact that the recombination process is non-planar. 

Based on these observations, one might expect that crosstalk processes and colour 
recombinations actually make very little numerical impact on the 2vl DPS cross section 
at small x. On the other hand, in an investigation of the four- gluon matrix element in the 
context of shadowing corrections to DIS [36], it has been shown that despite the 1/{N^ — 1) 
suppression of the 'recombination vertex', the inclusion of recombination effects in the 
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evolution of the four-gluon matrix element leads to an increase in this matrix element by 
approximately 70%, for values of x and an 'evolution length' t = Iii^Q^/Qq) relevant to the 
HERA experiment (x ~ 10~^, t ~ 3). According to equation (23) in [36] this correction 
should become even larger for smaller values of x and longer evolution lengths, such as 
those that will be relevant in DPS at the LHC. Therefore, we cannot simply ignore the 
effect of the recombination processes in the 2vl DPS cross section purely because of their 
colour suppression - they must be carefully taken into account. 

From the point of view of low x physics, there is an important distinction between 
the two crosstalk processes that we have discussed in this section - i.e. the off-diagonal 
real emission and virtual exchange processes. The off-diagonal real emission process can 
significantly reduce the magnitudes of the lightcone momentum fractions of the two active 
parton legs involved in the process, since it is a real emission process. On the other hand, 
the same is not true for the virtual exchange process. Here, the sum of the lightcone 
momentum fractions of the two parton legs involved must be conserved, and since the two 
legs are forced to have positive lightcone momentum fractions by the kinematics of the 
process, the magnitudes of both xs cannot simultaneously decrease - one must increase to 
compensate the decrease of the other. This means that, taking all partons involved to be 
gluons as is appropriate at low x, the virtual exchange splitting function is not enhanced 
at small x in the same way that the off-diagonal exchange (and indeed diagonal exchange) 
splitting functions are. In particular, the virtual exchange diagrams do not contribute at 
double leading logarithmic order to the evolution of the four-gluon matrix element. This 
result has been known for some time - see [33, 45-49]. This is the reason why we drew 
the colour recombination process in figure 6 using two off-diagonal real emission processes 
- it would be also possible to engineer a colour recombination using two virtual exchange 
processes instead, but such a process would not be as strongly enhanced at low x. 

The statement that the 'fo = 0' twist-4 distributions probed in the 2vl contribution 
to DPS evolve differently from the 2pGPDs with finite b has been made recently in Ap- 
pendix A of [24]. However, in this paper it is claimed that only the evolution of the colour 
correlation/interference distributions changes at 6 = - we contend that the evolution of 
the colour diagonal/singlet distribution is also affected in an important way. In equation 
(Al) of [24], an evolution equation for the colour octet twist-4 qq distribution diagonal in 
x fractions is proposed. However, the equation they propose involves only similar distribu- 
tions diagonal in x fractions on the right hand side - in fact the correct evolution equation 
should contain more general distributions nondiagonal in x on the right hand side, since the 
crosstalk processes that are allowed for 6 = will necessarily disrupt a diagonal/symmetric 
pattern of x values. 

4 The Total Cross Section for Double Parton Scattering 

In a previous publication [20] (see also [50]), we examined another class of diagrams that 
can contribute to the DPS cross section - the 'Ivl' graphs. For the arbitrary 'Ivl' graph in 
figure 1(c) with a total of n QCD branching vertices in the amplitude or conjugate, we dis- 
covered that there is no natural piece of the diagram that is proportional to log((5^/A^)"/i2p 
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and is associated with the transverse momenta inside the loop being strongly ordered on 
either side of the diagram (in fact, most of the contribution to the total cross section expres- 
sion for the graph comes from the region of integration in which the transverse momenta 
of particles inside the loop are of 0{\/Q^)). Based on this finding, we suggested that Ivl 
graphs should not contribute to the pp DPS cross section. 

Combining this suggestion with the findings of sections 2 and 3, we obtain the following 
formula for the total DPS cross section: 

D I \ D,2v2i \ , D,2vl, x 

'^{A,B)iV = («) + («) (4-1) 

with '^(^^^^(s) and o'^jf^^{s) being given by the expressions (2.34) and (3-14) respectively^. 

This formula agrees with the DPS cross section formulae proposed in [24, 26] and [19], 
apart from the fact that in neither of these papers are the crosstalk effects in the 2vl graphs 
taken into account correctly (in [19] they are omitted, whilst in [24, 26] they are included 
in an incorrect fashion). Ryskin and Snigirev [25] include an extra 'Ivl' contribution in 
their proposed cross section formula, which they argue can be validly included in the DPS 
cross section at suitably low x [32] , but which should not be included at moderate to large 
X values. 

We would like to point out at this stage that there are two features in the equation 
(4.1) that are potentially concerning, and that might indicate that modifications to it may 
be required in order to correctly describe the DPS cross section. 

The first issue is that we were originally expecting to obtain an expression for the DPS 
cross section looking something like (1.1), with the 2pGPDs in these formulae each having 
an interpretation in terms of hadronic operator matrix elements. Our proposed expression 
(4.1) deviates somewhat in structure from these expectations. 

The second issue is that there is a rather sharp distinction in (4.1) between perturba- 
tively and nonperturbatively generated parton pairs, with the 2pGPD for the latter having 
a natural width in r space of order A (as was discussed in section 2). Does there exist 
some scale at which we can (approximately) regard all parton pairs in the proton as being 
'nonperturbatively generated' in this sense (as is assumed in (4.1))? If so, what is the 
appropriate value for the scale (presumably it should be rather close to Aqc£))? 

These issues are related in an essential way to the fact that we have cut the contribution 
from 'Ivl' graphs out of the DPS cross section entirely. It may therefore not be correct to 
entirely remove these graphs from the DPS cross section in this way. On the other hand, 
at present we do not have a suitable alternative prescription for handling these graphs, and 
leave finding the appropriate way of including the Ivl graphs to future work. 

^Note that this is really our prediction for the unpolarised diagonal contribution to the total DPS cross 
section when the scales of the two hard interactions are the same, Q\ = Q% = . To generalise this result 
to unequal scales, one needs to change to Q\ in all Green's functions in (3.14) involving a '1' index, 
change to Q% in all Green's functions in (3.14) involving a '2' index, change the upper limit of the 

integration to mm{Q\,Q%), and perform a similar operation for the '2v2' contribution. As mentioned 
previously, the contributions associated with spin polarisation (either longitudinal or transverse) and flavour 
interference are expected to have a similar structure to (4.1), whilst the colour correlation/interference 
and parton type interference contributions should be suppressed by Sudakov factors, as is discussed in 
[5, 6, 26, 30]. 
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5 Conclusions 



In this paper we have closely examined the contribution to the LO p-p DPS cross section 
from graphs in which two 'nonperturbatively generated' ladders interact with two ladders 
that have been generated via a perturbative 1 — >• 2 branching process - '2vl' graphs. 
We have presented a detailed calculation demonstrating that 2vl graphs in which the two 
nonperturbatively generated ladders do not interact with one another contribute to the LO 
p-p DPS cross section in the way originally written down by Ryskin and Snigirev [25], and 
then later by Blok et al. [19] and Manohar and Waalewijn [24, 26]. We have also shown 
that 2vl graphs in which the 'nonperturbatively generated' ladders exchange partons with 
one another contribute to the LO p-p DPS cross section, provided that this 'crosstalk' 
occurs at a lower scale than the 1—7-2 branching on the other side of the graph. We have 
proposed a formula for the contribution from 2vl graphs to the LO DPS cross section, 
equation (3.14). 

Crosstalk interactions between the two nonperturbatively generated ladders are sup- 
pressed by colour effects, as has been noted in past studies of these interactions in the 
context of the twist 4 contribution to DIS. The most likely type of crosstalk interaction at 
low X is a 'colour recombination', in which the grouping of the four parton legs into two 
colour singlet ladders is altered - this is preferred because it maintains the overall system 
as two colour singlet ladders, and the anomalous dimension is largest for this configura- 
tion. For such a recombination, the colour factor associated with the vertex is 1/{N^ — 1). 
However, this suppression of the vertex does not necessarily mean that the effects of colour 
recombinations can be neglected - it was discovered in [36] that crosstalk effects signif- 
icantly affected the size of the four-gluon twist-4 operator for x values ~ 10~^ and an 
evolution length ln(Q^/Qo) — ^' ^'^ plausible that they could have an important im- 
pact in the context of DPS at the LHC. Clearly, more detailed investigations are required 
in order to assess the precise numerical contribution of colour recombinations on the DPS 
cross section. 

We combined our formula for the 2vl contribution to the DPS cross section (3.14) with 
the suggestion that we made in [20] that Ivl graphs should be completely removed from 
the DPS cross section to suggest a formula for the DPS cross section, equation (4.1). Two 
potentially concerning features were identified in this equation, and the existence of these 
might indicate that completely removing the Ivl graphs from the DPS cross section is not 
quite the correct prescription. The determination of the appropriate manner of treating 
the Ivl graphs is left to future work. 

Note Added 

After this paper was completed we learned of the published version of 'What is Double 
Parton Scattering?' by Manohar and Waalewijn [51] in which were corrected the errors of 
the arXiv version [24] that we discussed at the end of section 3. The discussion in that 
paper now appears to be in alignment with our own findings. 



-24- 



Acknowledgements 

We gratefully acknowledge discussions with J. Stirling, and financial support from a Trinity 
College Senior Rouse Ball Studentship. All the figures in this paper were generated using 
JaxoDraw [52]. 

References 

[1] A. Del Fabbro and D. Trcleani, A double parton scattering background to Higgs boson 
production at the LHC, Phys. Rev. D61 (2000) 077502, [hep-ph/99 11358]. 

[2] M. Y. Hussein, A double parton scattering background to associate WH and ZH production at 
the LHC, Nucl. Phys. Proc. Suppl. 174 (2007) 55-58, [hep-ph/0610207]. 

[3] M. Y. Hussein, Double parton scattering in associate Higgs boson production with bottom 
quarks at hadron colliders, arXiv : 0710 . 0203. 

[4] D. Bandurin, G. Golovanov, and N. Skachkov, Double parton interactions as a background to 
associated HW production at the Tevatron, JHEP 1104 (2011) 054, [arXiv : 1011 . 2186]. 

[5] M. Diehl and A. Schafcr, Theoretical considerations on multiparton interactions in QCD, 
Phys.Lett. B698 (2011) 389-402, [arXiv : 1102 . 3081]. 

[6] M. Diehl, D. Ostermeier, and A. Schafer, Elements of a theory for multiparton interactions 
in QCD, JHEP 1203 (2012) 089, [arXiv : 1111 . 0910]. 

[7] Axial Field Spectrometer Collaboration, T. Akesson et. al.. Double parton scattering in 
pp collisions at y/^ ^ 63GeV, Z. Phys. C34 (1987) 163. 

[8] UA2 Collaboration, J. Alitti et. al., A Study of multi - jet events at the CERN anti-p p 
collider and a search for double parton scattering, Phys.Lett. B268 (1991) 145-154. 

[9] CDF Collaboration, F. Abe et. al.. Study of four jet events and evidence for double parton 
interactions in pp collisions at ^/s = 1.8 TeV, Phys. Rev. D47 (1993) 4857-4871. 

[10] CDF Collaboration, F. Abe et. al.. Double parton scattering in pp collisions at 
yi= 1.876^^, Phys. Rev. D56 (1997) 3811-3832. 

[11] DO Collaboration, V. M. Abazov et. al.. Double parton interactions in photon+3 jet events 
in p p- bar collisions sqrts=L96 TeV, Phys. Rev. D81 (2010) 052012, [arXiv : 0912 . 5104]. 

[12] A measurement of hard double-partonic interactions inW-^lv + 2 jet events, Tech. Rep. 
ATLAS-CONF-2011-160, CERN, Geneva, Dec, 2011. 

[13] P. Bartalini (ed.) and L. Fano (ed.). Multiple partonic interactions at the LHC. Proceedings, 
1st Lnternational Workshop, MPL'08, Perugia, Italy, October 27-31, 2008, arXiv : 1003 . 4220. 

[14] http://nipilO.desy.de/. 

[15] http://www.mpi2010.physics.gla.ac.uk/Home.htnil. 
[16] http://nipill.desy.de/. 

[17] P. Bartalini, E. Berger, B. Blok, G. Calucci, R. Corke, et. al., Multi-Parton Interactions at 
the LHC, arXiv: 1111.0469. 

[18] B. Blok, Y. Dokshitzer, L. Frankfurt, and M. Strikman, The Four jet production at LHC and 
Tevatron in QCD, Phys. Rev. D83 (2011) 071501, [arXiv : 1009 . 2714]. 



-25- 



[19] B. Blok, Y. Dokshitzer, L. Frankfurt, and M. Strikman, pQCD physics of multiparton 
interactions, Eur.Phys.J. C72 (2012) 1963, [arXiv : 1106 . 5533]. 

[20] J. R. Gaunt and W. J. Stirling, Double Parton Scattering Singularity in One- Loop Integrals, 
JHEP 1106 (2011) 048, [arXiv : 1103 . 1888]. 

[21] V. P. Shclcst, A. M. Snigirev, and G. M. Zinovjcv, Gazing into the multiparton distribution 
equations in QCD, Phys. Lett. B113 (1982) 325. 

[22] G. M. Zinovev, A. M. Snigirev, and V. P. Shelest, Equations for many parton distributions in 
quantum chromodynamics, Theor. Math. Phys. 51 (1982) 523-528. 

[23] A. M. Snigirev, Double parton distributions in the leading logarithm approximation of 
perturbative QCD, Phys. Rev. D68 (2003) 114012, [hep-ph/0304172]. 

[24] A. V. Manohar and W. J. Waalewijn, What is Double Parton Scattering?, arXiv : 1202 . 5034. 

[25] M. Ryskin and A. Snigirev, A Fresh look at double parton scattering, Phys. Rev. D83 (2011) 
114047, [arXiv: 1103. 3495]. 

[26] A. V. Manohar and W. J. Waalewijn, A QCD Analysis of Double Parton Scattering: Color 
Correlations, Interference Effects and Evolution, arXiv : 1202 . 3794. 24 pages, 15 figures. 

[27] N. Paver and D. Treleani, Multi-quark scattering and large px jet production in hadronic 
collisions, Nuovo Cim. A70 (1982) 215. 

[28] M. Mekhfi, Multiparton processes: an application to double Drell- Yan, Phys. Rev. D32 
(1985) 2371. 

[29] M. E. Peskin and D. V. Schroeder, An Introduction to quantum field theory. Addison- Wesley 
Publishing Co., 1995. 

[30] M. Mekhfi and X. Artru, Sudakov suppression of color correlations in multiparton scattering, 
Phys. Rev. D37 (1988) 2618-2622. 

[31] J. R. Gaunt and W. J. Stirling, Double Parton Distributions Incorporating Perturbative QCD 
Evolution and Momentum and Quark Number Sum Rules, JHEP 1003 (2010) 005, 
[0910.4347]. 

[32] M. Ryskin and A. Snigirev, Double parton scattering in double logarithm approximation of 
perturbative QCD, arXiv: 1203.2330. 

[33] E. Levin, M. Ryskin, and A. Shuvaev, Anomalous dimension of the twist four gluon operator 
and pomeron cuts in deep inelastic scattering, Nucl.Phys. B387 (1992) 589-616. 

[34] J. Bartels, Unitarity corrections to the Lipatov pomeron and the small x region in deep 
inelastic scattering in QCD, Phys. Lett. B298 (1993) 204-210. 

[35] J. Bartels, Unitarity corrections to the Lipatov pomeron and the four gluon operator in deep 
inelastic scattering in QCD, Z.Phys. C60 (1993) 471-488. 

[36] J. Bartels and M. Ryskin, Absorptive corrections to structure functions at small x, Z.Phys. 
C60 (1993) 751-756. 

[37] J. Bartels and M. Ryskin, Recombination within multi-chain contributions in pp scattering, 
arXiv: 1105. 1638. 

[38] R. K. Ellis, W. Furmanski, and R. Petronzio, Power Corrections to the Parton Model in 
QCD, Nucl.Phys. B207 (1982) 1. 

[39] R. K. Ellis, W. Furmanski, and R. Petronzio, Unraveling Higher Twists, Nucl.Phys. B212 



-26- 



(1983) 29. 

[40] J.-W. Qiu, Twist four contributions to the parton structure functions, Phys.Rev. D42 (1990) 
30-44. 

[41] J.-W. Qiu and G. F. Stcrman, Power corrections in hadronic scattering. 1. Leading \/Q^ 
corrections to the Drell-Yan cross-section, Nucl.Phys. B353 (1991) 105-136. 

[42] A. Bukhvostov, G. Frolov, L. Lipatov, and E. Kuraev, Evolution Equations for 
Quasi- Partonic Operators, Nucl.Phys. B258 (1985) 601-646. 

[43] L. Gribov and M. Ryskin. LNPI-865 (1983). 

[44] L. Gribov, Y. Dokshitzcr, S. Troyan, and V. Khozc Sov. Phys. JETP 88 (1988) 1303. 

[45] B. M. McCoy and T. T. Wu, Mandelstam Diagrams Are Not Enough, Phys.Rev. D12 (1975) 
546. 

[46] B. M. McCoy and T. T. Wu, Three Particle Regge Pole in (f>^ Theory, Phys.Rev. D12 (1975) 
578. 

[47] S. Matinyan and A. Scdrakyan JETP Lett. 23 (1976) 588. 
[48] S. Matinyan and A. Scdrakyan JETP Lett. 24 (1976) 240. 
[49] S. Matinyan and A. Scdrakyan 5*0?;. J. Nucl. Phys. 24 (1976) 844. 

[50] J. R. Gaunt and W. J. Stirling, Single and Double Perturbative Splitting Diagrams in Double 
Parton Scattering, arXiv: 1202 . 3056. 

[51] A. V. Manohar and W. J. Waalcwijn, What is Double Parton Scattering?, Phys. Lett. B713 
(2012) 196-201. 

[52] D. Binosi and L. Thcufil, JaxoDraw: A graphical user interface for drawing Feynman 
diagrams. Computer Physics Communications 161 (2004) 76-86. 



-27- 



